Flat-lens focusing of electrons on the surface of a topological insulator 
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We propose the implementation of an electronic Veselago lens on the conducting surface of a three- 
dimensional topological insulator (such as Bi2Te3). The negative refraction needed for such a flat 
lens results from the sign change of the curvature of the Fermi surface, changing from a circular to a 
snowfiake-like shape across a sufficiently large electrostatic potential step. No interband transition 
(as in graphene) is needed. For this reason, and because the topological insulator provides protection 
against backscattering, the potential step is able to focus a broad range of incident angles. We 
calculate the quantum interference pattern produced by a point source, generalizing the analogous 
optical calculation to include the effect of a noncircular Fermi surface (having a nonzero conic 
constant). 
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I. INTRODUCTION 

Ballistic electron optics relies on the analogy be- 
tween the Schrodinger equation for electrons and the 
Helmholtz equation for classical waves to construct de- 
vices that can image the flow of electrons in high-mobility 
semiconductors. 1_4 A variation in electrostatic poten- 
tial is analogous to a variation in dielectric constant, 
so that a curved gate electrode can have the refractive 
power of a lens — as has been demonstrated in the two- 
dimensional electron gas of a GaAs heterostructure. 5,6 
The focal length of this electrostatic lens depends on its 
curvature, diverging for a flat electrode. 

Focusing of light by a flat lens is possible in media 
with a negative index of refraction. This so-called Vese- 
lago lens 7 ' 8 has a focal length proportional to the distance 
between lens and source, rather than fixed by the lens it- 
self. It is also not limited by the single optical axis of a 
curved lens and can have a much wider aperture. Pho- 
tonic crystals can provide the negative refraction needed 
for a flat lens, 9 as demonstrated experimentally. 10,11 

The electronic analog of a Veselago lens was proposed 
in the context of graphene, 12 based on the negative re- 
fraction of an electron crossing from the conduction band 
into the valence band. Such interband crossing requires 
a p-n junction, which is highly resistive if the interface 
extends over more than an electron wave length. 13,14 It 
would be desirable to have a method for producing a 
flat lens entirely within the conduction band, in order to 
avoid a resistive interface. It is the purpose of this work 
to propose such a method, in the context of topological 
insulators. 

Topological insulators have a conducting surface with 
a Dirac cone of massless, helical low-energy excitations, 
reminiscent of graphene. 15 ' 16 Indeed, scanning tunneling 
microscopy has shown that backscattering of the surface 
electrons is inhibited, as expected from conservation of 
helicity. 17 ' 18 While the large band gap topological insula- 
tor Bi2Se3 has a nearly circular Dirac cone, in the smaller 
band gap material Bi2Te3 the cone is warped in an hexag- 
onal snowflake-like shape. 19 ~ 23 The hexagonal warping 



of the Fermi surface enhances the quantum interference 
(Friedel) oscillations in the density of states near an im- 
purity or potential step, 24-27 which for a circular Fermi 
surface would be suppressed by conservation of helicity. 28 

The electron focusing considered here is an altogether 
different, semiclassical consequence of the hexagonal 
warping. The flat lens is formed by a potential step on 
the surface of the topological insulator, sufficiently high 
to change the curvature of the Fermi surface from con- 
vex to concave. The sign change of the curvature leads to 
negative refraction and focusing, qualitatively similar to 
the optical Veselago lens — but quantitatively different 
because of the nonuniformity of the curvature (quantified 
by a nonzero conic constant). 

In the following two sections, we derive the negative 
refraction and the line of focal points (caustics), as well 
as the diffraction pattern produced by a point source. We 
calculate the curvature and conic constant for the specific 
case of Bi2Te3. We conclude in Sec. IV by comparing 
with the flat lens formed by a p-n junction in graphene 12 
and by discussing possible experimental realizations in 
topological insulators. 



II. NEGATIVE REFRACTION AT A 
POTENTIAL STEP 

A. Negative refraction 

Consider an electron propagating approximately along 
the x-axis (the optical axis) and impinging at x = onto 
an electrostatic potential step 5U produced by a gate 
electrode (see Fig. 1). For simplicity, we assume that the 
optical axis is parallel to an axis of crystallographic sym- 
metry, such that the equi-energy contours are ±k y sym- 
metric. (For the more general case, see App. A.) At con- 
stant Fermi energy, the kinetic energy changes from Ei in 
the incident (left) region to E t = Ei + SU in the transmit- 
ted (right) region. The equi-energy contour at the left is 
given locally by 5ki. x = — \c,ik\ y , for a two-dimensional 
wave vector ki = {k%fi + 8ki tX , ki. y ) approximately along 




FIG. 1: Negative refraction at a potential step (height SU) 
where the curvature of the equi-energy contours (thin curves) 
changes sign from a > to c t < 0. An electron (thick arrow) 
with kinetic energy Ei is incident at angle 6i and transmitted 
at angle 6t- Because the curvature changes sign, the electron 
is negatively refracted with 9 t < for 0; > 0. 



the optical axis, and similarly Skt jX = — -^Ctkfy at the 
right. The coefficients Cj and c t are the curvatures of the 
Fermi surface for normal incidence, at the two sides of 
the potential step. 

The velocity v = ftr 1 dE/dk is normal to the equi- 
energy contours, so that the velocities Vi and Vt in the 
left and right regions make, respectively, an angle 9i = 
Ciki^y and 9 t — Ctk t ^ v with the x-axis. Conservation of 

k t .y) leads to the linearized 



transverse momentum (ki^ y 
Snell's law 

6 t = (ct/cOft, 



for 6 u 6 t < 1. 



(2.1) 



The inverse curvature plays the role of the refractive in- 
dex in optics. Negative refraction (meaning 6i$t < 0) 
takes place when d and ct have opposite signs, as illus- 
trated in Fig. 1. 



B. Noncircular Snell's law 

As we will see in the next section, to calculate the 
image of a point source we will need to include the first 
nonlinear correction to Eq. (2.1). In optics, where one 
has a circular equi-energy contour, Snell's law Cj sin# t = 
ct sin 9i implies the series expansion 



with m — Ct/ci and n 3 
we can write 



e t = nidi + n 3 ef + o{e!), (2.2) 

hni{n\ — 1). More generally, 



n\ = Ct/cu n 3 = \ni{n\ - 1) + A, 



(2.3) 



where A quantifies the deviation from the optical Snell's 
law. 31 

The parameter A vanishes for a circular Fermi sur- 
face as in graphene, 12 ' 29,30 but is nonzero for the warped 
Fermi surfaces of topological insulators. In order to relate 
A to the Fermi surface, we parameterize the equi-energy 



contour using polar coordinates by k = k(</>)(cos0, sin0), 
where (j> is the angle between the wave vector k and the 
x-axis and k = \k\. A subscript i or t distinguishes the 
parameters at the two sides of the potential step. 

The noncircular Snell's law is expressed by the three 
equations 



K t (4> t ) sin (f> t = Ki(4>i) sin fa , 

K, (</>;) tan 0; - Ki((j)i 



tanw, = 



tan#t 



K t ((f>t)ta,ii(j>t - K' t ((f> t ) 
K t (<f>t) + «J((/) t )tan0 t ' 



(2.4) 
(2.5) 

(2.6) 



where k' — dn/dcj). The first equation expresses the con- 
tinuity of the y-component of the wave vector at the in- 
terface x = 0, while the second and third equations relate 
the angles 9 and <j> of velocity and wave vector (using the 
fact that v is perpendicular to the equi-energy contour). 
The circular Snell's law K t sm9 t = K^sin^ is recovered 
for k' — 0, when 9 — (p. 

Near <f> = 0, the equi-energy contour can be param- 
eterized in terms of the curvature c and conic constant 
K, 

k 2 y = -(2/c)Sk x -(l+IC)(Sk x ) 2 , (2.7) 

with 8k x — k x — k(0). The noncircular Snell's law then 
expands to Eqs. (2.2)-(2.3) with 



Q 



* - £(&. 



CjICi 



(2.8) 



C. Application to Bi2Te3 

We apply these general considerations to the topolog- 
ical insulator Bi2Tc3. On the [111] surface and close to 
the center of the Brillouin zone (the T point) the Hamil- 
tonian can be approximated by 28 

H = Hv k (ay cos 4> — o x sin <p + A k a z cos3<fi) . (2.9) 

The dispersion relation in the conduction band (E > 0) 



E(k) = Wfc 2 + (A 2 fc 3 cos3</>) 5 



= hvy / k 2 + \ 4 (k 3 x -3k x kl) 



2^2 
V> 



(2.10) 



The <Tj's are Pauli matrices acting on the electron spin 
and 4> denotes the angle of the wave vector k with respect 
to the TK direction in the Brillouin zone (oriented along 
the x-axis). The parameters v « 4 • 10 5 m/s and A rj 
1 nm were estimated by fitting to data from angularly 
resolved spectroscopy. 21,22 (Additional terms quadratic 
in momentum can be included in the fit, but these do 
not qualitatively change the dispersion.) 

The curvature c(E) of the equi-energy contour in the 
TK direction is given by 



c(E) = A 



l-6^ 4 
x + 3>r 5 ' 



(2.11) 
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FIG. 2: Curvature c(E) of the equi-energy contour in the 
TK direction, calculated from Eq. (2.11). The shape changes 
from convex to concave at energy E c . The maximally negative 
curvature is c ~ —1.3 A for E ~ 2 hv/X, where the equi-energy 
contour has the snowflake-like shape shown in the inset. 




EX/hv 



FIG. 3: Plot of the conic constant /C (solid line) as well as 
the combination c 3 /C/A 3 [appearing in the noncircular Snell's 
law (2.8)] (dashed line), both as a function of the energy E. 
The divergence of K, is at the energy E c where the curvature 
vanishes. 



with x 2 = £ , — £_ defined in terms of 



e± = 



6h 2 v 2 



§?i 4 i; 4 + 9A 4 £ 4 ± 3A 2 £ 2 ) . (2.12) 



The quantity x/A = /c(0) equals |fe| at <j> = 0. 

The energy dependence of the curvature is plotted in 
Fig. 2. As discovered by Fu, 28 the curvature changes 
sign when x* — 1/6, which corresponds to an energy 
E c = 6~ 3 / 4 \/7 hv/X ss 0.2 eV and a wave vector k c = 
x c /X sa 0.6nm _1 . At the same point the conic constant 



JC(E) = 



3x 4 (35-60x 4 + 72x 8 



(1 - 6x 4 ) 3 
diverges and thereby changes sign, cf. Fig. 3. 



(2.13) 




FIG. 4: Classical trajectories refracted at a potential step at 
x — 0, with the cusp caustic indicated. 



III. CAUSTICS FROM A POINT SOURCE 
A. Focusing of classical trajectories 

Because of the negative refraction, diverging trajecto- 
ries become converging at the potential step and then 
cross at a focal point (see Fig. 4). If a point source 
is placed at (— a, 0), a distance a from the interface at 
x = 0, then the trajectory for an electron incident at an 
angle 9i and transmitted at an angle 8 t is parameterized 

by 



y{x;8i) 



(a + x) tan#i, 
atan(9j + x tan#t, 



for x < 0, 
for x > 0. 



(3-1) 



On the optical axis y, 
(of,0), with 



,^->0we obtain the focal point 



Of 



-a/ni 



— a, 

Ct 



(3.2) 



proportional to the ratio of the two curvatures. As in 
the optical Veselago lens, 32 ' 33 the focal point is displaced 
from the optical axis as we increase the angle of incidence, 
so that the point (aj?, 0) is the cusp on a curve of focal 
points. This caustic curve (called an astroid 34 ) is visible 
in Fig. 4 as the envelope of the refracted trajectories. 

The caustic curve near (aF,0) is obtained from Eq. 
(3.1) and the nonlinear Snell's law (2.2) by demanding 
that dy/dOi = 0. We find 



a{y/a) 2 = {x/a F - if 



(3.3) 



with the opening rate of the cusp governed by the pa- 
rameter 



27 



a = — [1 - {ct/ci) 2 - 2(ci/ct)A] 



(3.4) 



For A = 0, so for a circular Fermi surface, this agrees 
with Refs. 12,33. Depending on the sign of a, the cusp 
points away from the potential step (for a > 0) or to- 
wards the potential step (for a < 0). For a — higher 
than third-order terms in the expansion (2.2) have to be 
included in order to obtain the caustic curve. 
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FIG. 5: Gray-scale plot of the current density j(r) as a func- 
tion of position r = (x,y) near the focal point (a,F,0) for 
a source located at (— a, 0), calculated from Eqs. (3.7) and 
(3.8) for a — 1 and a/ci — 100. Lighter shades of gray in- 
dicate higher current densities. The cusp caustic starting at 
the focal point is decorated by oscillations on the scale of the 
wavelength. 



B. Quantum interference near the focal point 



The diffraction pattern near a cusp caustic has a uni- 
versal functional form (Pearcey integral), 35 ' 36 but the pa- 
rameters governing that function are modified for non- 
circular equi-energy contours. We calculate the wave 
function f at a point r = (x, y) near the cusp by sum- 
ming over partial waves ^ yo from points r$ — (0,t/o) 
along the potential step [excited by a point source at 
^source = (—0,0)]. In the far-field approximation, for a 
and cif large compared to the wave length, the partial 
waves have the simple form 



*!/«,= 



A e^y 



$ Vo = ki ■ (r - r s 



k t -{r- r ) 



(3.5) 
(3.6) 



The amplitude A ya and spinor components u yo , v yo vary 
slowly as yo is varied on the scale of the wave length, so 
we fix their values at Ao,uo,vo and retain only the yo 
dependence of the phase <& ya . 

In the optical case, the wave vectors ki and k t at the 
two sides of the interface point in the direction of the 
velocity and hence are parallel to the rays r$ — r source 
and r — rr,. For a noncircular Fermi surface this is no 
longer true and we have to take into account the differ- 
ence between the angles cf>i, 4>t and 9i = arctan(yo/a) 
and 9 t — arctan [(y — yo)/x] which the wave vectors and 
the rays make with the x-axis. The relation between <f> 
and 9 is expressed by Eqs. (2.5) and (2.6), in terms of the 
radial parameter k(4>) = \k\ of the equi-energy contour. 

We expand <& yo in a power series in yo- Near the 
cusp caustic (3.3), y/a = O(yo/a) 3 while x/ap — 1 = 




FIG. 6: Flat lens, with two potential steps (upward at x = 
and downward at x = L) and two cusp caustics (at x — a\ 
and x — ap). 



0(yo/ciF) 2 . To fourth order in yo we find 



%o = n t (Q)a + K t (0)x 

yyo 



(x - a F )y 2 



aj/o 



CtCiF ^cto-p 27 Cta fo? 



(3.7) 



with a given by Eq. (3.4). One readily checks that the 
stationary phase equations d$ yo /dyo = = d 2 & Vo /dy 1 
give the caustic curve (3.3). (These equations correspond 
to the geometric optics limit c t — > of vanishing wave- 
length.) 

The current density j (r) follows upon integration over 

yo, 

2 



j( r ) 



.70 



dy e 



'*./ 



(3.8) 



with jo a constant proportional to the product of the 
injection rate at the source and the transmission prob- 
ability T through the potential step. By rescaling the 
integration variable y — > ay , we see that the current 
density (3.8) as a function of x/ap and y/a depends only 
on the two parameters a and a/ci. Figure 5 is a plot of 
this current density, showing the characteristic interfer- 
ence pattern of a cusp caustic. 



C. Focusing by a fiat lens 

The flat lens in Fig. 6 is formed by the potential profile 
U(x) = SU for < x < L, U(x) = otherwise. We 
denote the Fermi surface curvatures (of opposite sign) 
inside the lens (0 < x < L) by ci cns and outside [x < 0, 
x > L) by co- Negative refraction at the two potential 
steps at x = and x = L focuses a source at x = —a 
on the optical axis (y = 0) first onto the point ai = 
— (co/ci cns )a inside the lens and then onto the point 

a F = (1 - ci cns /c )L - a (3.9) 

outside the lens (provided it is sufficiently thick, 
|ciens£| > |c a|). 
The classical trajectories are now parameterized by 

{(a + x)ta,n6i, for x < 0, 
a tan 9i + x tan 9 t , for < x < L, 
(a + x — L) tan 9i + L tan. Ot, for x > L. 



(3.10) 



The relation between 9 t and 9% is still given by Eq. (2.2), 
with 



n\ = c lons /c , n 3 = \n\{n\ - 1) + A. 
The cusp caustic near (ap-,0) has the form 

(i(y/a) 2 = (x/a F - if, 
as in Eq. (3.3) but with a different parameter 

27 La 2 cie 



(3.11) 



(3.12) 



P 



CO 



1 - (cie„ s /co)' - 2(c /c lens )A . (3.13) 



Notice that a and j3 have the opposite sign (because of 
the factor c\ cns /cQ < 0), so that the cusps inside and 
outside the lens point in opposite directions (as visible in 
Fig. 6). 

The flat-lens diffraction pattern near the caustic is 
given by the same Pearcey integral (3.7)-(3.8) as for a 
single interface, but with different coefficients, 



j{r) = jo 



dy exp 



yyo 

c a F 



(x - a F )y 2 
2c a| 



27 cqQf a, 2 



(3.14) 



Thus, the interference pattern that can be observed near 
a F looks similar to Fig. 5. 



IV. DISCUSSION 

A. Intraband versus interband negative refraction 

The Veselago lens at a p-n junction in graphene 12 uses 
interband scattering to achieve negative refraction. In 
contrast, the mechanism considered here is intraband, 
operating entirely within the conduction band. The p- 
n junction has one special feature which our setup lacks, 
which is the possibility to use electron-hole symmetry 
to collapse the caustic curve onto a single focal point 
(when c t = —Ci). In our setup the Fermi surfaces at the 
two sides of the potential step are not related by any 
symmetry relation, so in general the two Fermi surface 
curvatures c t and a will be different in magnitude. 

The main advantage of an intraband over an inter- 
band mechanism for negative refraction is that the trans- 
mission probability T can be much higher. Typically, 
the width d of the potential step will be large com- 
pared to the Fermi wave length X F = 2 r K/k F . Intra- 
band transmission is then realized with unit probabil- 
ity, up to exponentially small backscattering corrections: 
T = 1 — 0(e^ kpd ). Interband transmission, in con- 
trast, has T ~ exp(— fc^iisin 9i), so it is exponentially 
suppressed for angles further than ^JXp/d from normal 
incidence. 14 



B. Experimental realization 

Realization of the intraband flat lens proposed here, 
requires firstly a topological insulator with sufficiently 
long mean free paths to ensure ballistic motion of the 
electrons from source to focus. Sufficiently pure single 
crystals should make this possible. 

Secondly, and more specifically, the curvature of the 
Fermi surface should be tunable from positive to nega- 
tive values by a gate voltage. From spectroscopic data 22 
for Sn-dopcd Bi2Te3 we would estimate that a potential 
step SU ~ — 0.1 eV would produce a positive curvature 
inside a narrow strip and a negative curvature outside 
(as in Fig. 6). The strip itself would also allow for bulk 
conduction, because in Bi2Te3 a positively curved Dirac 
cone of surface states overlaps with bulk states. Since the 
regions outside the lens have only surface conduction, we 
do not expect the bulk states inside the lens to spoil the 
focusing. 

A point source can be created, for example, using 
the "needle-anvil" technique developed for point contact 
spectroscopy, 37 or alternatively using a scanning tunnel- 
ing microscope (STM). For the spatially resolved detec- 
tion of the current density distribution an STM tip is 
most convenient. Such a setup would provide a sensitive 
probe of the nonspherical Fermi surface of a topological 
insulator, in a similar way as has recently been proposed 
for metals. 38 
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Appendix A: Sheared caustic curve for tilted 
potential step 

In the main text we have assumed for simplicity that 
the potential step is perpendicular to the TK direction 
in Fig. 2. Then only odd powers of 9i appear in the 
expansion (2.2). If the potential step is tilted relative to 
this crystallographic axis, then the cusp caustic persists 
but in a distorted form, as we now derive. 

Including also even powers of 9i in Eq. (2.2) one would 
have the expansion 



l 2 u i 



l 3 u i 



0(9j 



(Al) 



By rotating the coordinate axis, we can set n' — 0. The 
expressions simplify if we expand in powers of tan^, 

tan 9t = mi tan 9% + m-z tan 2 6i + m 3 tan 3 6i + C(tan 4 6i). 

(A2) 



From Eq. (3.1), demanding dy/dOi = 0, we obtain the 
implicit caustic equation 



a(mi + 2?ri2 tan Bi + Sm^ tan 6i) 



-1 
(777,2 + 2rn3 tan $i) tan 2 (9, 



(A3) 



The cusp of the caustic is given by the condition 
dx/dOi = 0. It is at tan# i0 = — to 2 /3to3. In order to 
remain in the region of validity of the expansion (Al), 
we assume that \m2\ <C |tos| so that the tilt remains 
small. Then the cusp is located at 



-a/mi 




(A4) 



We now expand Eq. (A3) near tan^ = tan^o to third 
order in 5 = tan 9i — tan 8io , 



X - x a \ _ | 

y-yo ' 



3mi<5 2 
-mim~l 1 S 2 + 2d 3 



(A5) 



Eliminating S yields the caustic curve 



i[y -Vo + e(x - x )} 2 = (x - x ) 3 , 



(A6) 



with coefficients 7 = 21am^/Am\ and e = 7771777,2/3777.3. 
Equation (A6) has the general form of a sheared cusp 
caustic from catastrophe theory. 39 
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